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Abstract 

^ \ The proton-neutron monopole Lipkin model, which exhibits some 

properties which are relevant for those double beta decay {(5(3) tran- 
sitions mediated by the Fermi matrix elements, is solved exactly in 
^ ', the proton-neutron two-quasiparticle space. The exact results are 

^ ' compared with the ones obtained by using the Quasiparticle Ran- 

dom Phase (QRPA) and renormalized QRPA (RQRPA) approaches. 
It is shown that the RQRPA violates the Ikeda Sum Rule and that 
this violation may be common to any extension of the QRPA where 
scattering terms are neglected in the participant one-body operators 
as well as in the Hamiltonian. This finding underlines the need of 
additional developments before the RQRPA could be adopted as a 
reliable tool to compute (3(3 processes. 
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1 Introduction 



The nuclear double beta decay could provide evidence on the existence of 
massive neutrinos and right handed weak currents [|IJ. This exciting possi- 
bility has attracted many theoretical and experimental work in the last years 
||. In order to extract information about this new physics the data must be 
complemented with theoretical nuclear matrix elements which are strongly 
suppressed. 

Ten years ago it was realized that the quasiparticle random phase approx- 
imation (QRPA), including a particle-particle channel in the residual interac- 
tion, can reproduce the experimentally determined two-neutrino double-beta 
decay (/3/3 2 ^) half-lives ||, |], |J, ||. This step allowed the description of sin- 
gle and double-beta decay transitions in many nuclei. However it was soon 
recognized that the ground state to ground state fifav transition amplitudes 
are extremely sensitive to the force parameters, thus limiting the predictive 
power of the theory. The breakdown of the QRPA approach, for some criti- 
cal values of the model parameters, made the theoretical description of some 
cases particularly difficult; i.e. the (3f3 2v decay of 100 Mo §, |, 0, g, §. 

Recently the use of a correlated vacuum in the QRPA equation of motion, 
the so-called renormalized QRPA (RQRPA) jrUL O, has been reformulated 



[ |12jl and applied to the Pfov decay |13|]. It was found that the formalism is 
stable beyond the point where the QRPA collapses. The RQRPA method 
requires the solution of coupled-non-linear equations, instead of the usual 
eigenvalue problem of the QRPA. However the physical consequences of this 
highly non-linear behaviour, represented by the inclusion of some terms be- 
yond the QRPA order of approximation, have not been explored carefully. 
In the present paper we will use a simple solvable model, which is an 



extension of the Lipkin model fll4| , to compare the exact, the QRPA and the 
RQRPA approaches. It will be shown that if one remains at the leading order 
of approximation, i.e. by retaining two-quasiparticle terms in the relevant 
Fermi or Gamow- Teller transition operators as well as in the Hamiltonian, 
the RQRPA violates the Ikeda sum rule. 
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2 The model 



The model hamiltonian |15], [16], [17j consists in a single particle, a pairing 
term for protons and neutrons and a schematic charge- dependent residual 
interaction including particle-hole and particle-particle channels. It has been 
shown in a recent series of papers that this interaction, treated in the frame- 
work of the QRPA, is as good as a G-matrix constructed from the OBEP 
Bonn potential in reproducing single- and double-beta decay matrix elements. 

The schematic hamiltonian reads 

H = H p + H n + H res , (1) 

where 

Hp = J2 e p a£a p — GpS^Sp H n = E e n t4a n — G n SlS n 

p n (2) 

H res = 2 X (3J ■ (3j - 2kPJ ■ Pj 

In the above expression the following definitions were introduced 



p n 

P7-ft= E (-1) M : PJmWJ-mV ■■ 

M=-J 

Pj-Pt= E (-1) M : PJm(PJ-mV ■ ( ' J>) 

M=-J 

Pjm = Ey < A°jm\j > ajdj PJ M = Eij < i\Oj M \j > a\aX 
Oim = vmt~ Coo = t~ 

a p = Oj mp being the particle creation operator and ap- = (—iy p ~ mp aj p _ m its 
time reversal. 

As mentioned above, the hamiltonian ([I]) with J = 1 provides a reason- 
able description of the main physics involved in Gamow- Teller transitions. 
The parameters x an d k play the role of the renormalization factors g p h and 
g pp introduced in the literature || f|, ||. § . 

For the case of J = the hamiltonian ([!]) can be reduced to an isospin 
scalar if its parameters are selected as 
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e p = e n , x = 0, G p = G n = 4k. (4) 

If x 7^ the isospin symmetry is broken in the particle-hole channel and 
for 4k 7^ G p , G n this symmetry is broken in the particle-particle channel. A 
similar identification of this isospin breaking mechanism can be performed 
for the case of realistic interactions with renormalized proton-proton (gp air ), 
neutron-neutron (g pair ) and proton-neutron particle-particle (g pp ) strengths 
[H, |4|, H P). For the Fermi transitions, like the ones contributing to the 
nuclear matrix elements associated with the neutrinoless double-bet a- decay 
mode, the value g pp = 1.0 is usually adopted, while g pair and g^ aiT vary from 
0.9 to 1.2 to reproduce the observed proton and neutron pairing gaps in 
medium and heavy mass nuclei. This parametrization, adopted for realistic 
interactions, provides a first motivation for using the proton-neutron particle- 
particle strength parameter independent of G. 

In order to construct an exactly solvable model to be used to test the reli- 
ability of the RQRPA solutions for the single- and double-beta decay observ- 
ables, we will make a strong approximation. We will consider a single shell 
with the same angular momentum for protons and neutrons, i.e. j p = j n = j. 
Adopting this model space is equivalent to work with the one-level limit of 
([]]). It will be shown that this model, which is not intended to reproduce 
actual nuclear properties, does have the qualitative features of a realistic 
pn-QRPA calculation. Indeed, excitation energies, single- and double-beta 
decay transition amplitudes and ground state correlations depend on the 
particle-particle strength parameter k in the same way as they do in more 
elaborated calculations with many single particle levels and with more real- 
istic interactions. The advantage of such a simplification lies upon the fact 
that the proton-neutron excitations can be described within the framework 
of an exactly solvable model. Particularly, the correspondence between the 
simplified version (one-shell limit) of (|I|) and the ordinary monopole Lipkin 
model |14j] can be established if one chooses the channel with J = of hamil- 



tonian Physically, this case will correspond to Fermi type transitions 
but we should emphasize the fact that the study of the model and not the 
adjustment of a given decay channel constitutes the essential aspect of the 
present discussion. 

In the single shell case, and after performing the Bogolyubov transforma- 
tions, separately, for protons and neutrons, the first two terms in ([T]) become 
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diagonal, i.e. 



Hp G p a pm p a pm p i H n 



(5) 



being e p , e n the quasiparticle energies and a P , c>4 the quasiparticle creation 
operators |20 



The linearized hamiltonian, neglecting the scattering terms (a p a ni a^ap) 
which do not contribute at the QRPA order, reads 



H € p a pm p a pm p f e nJ2 a nm,r, a 



+ 



2 X (2j + 1) [( U y n + vyjA ■ a + 

UpV n v p u n A^ ■ + VpU n u p v n A ■ A 



(6) 



2«(2j + l) [(uy n + vy n )A^A- 

UpU n v p v n A^ ■ A^ — VpV n u p u n A • A 



with 



At 



j=o 

M=0 



To show that it is possible to reduce (||) to an holomorphic version of the 
Lipkin model we have considered, for simplicity, one single particle orbital 
for protons and neutrons with the same quasiparticle energies e p = e n = 
(2j + 1)^/4, j = P) n0. 

Under these approximations the hamiltonian (0) has the form 



H = eC + XiA^A + \ 2 {A ] A ] + AA) 



(7) 



with 



C — S a pm v a pm p + a nm„ a nm l ■ 
rn p r m n 



and 



Ai 



X(uy + vy n ) - n(u 2 y n + vy n ) , A 2 = 2(x + K)u p v p u n v 



jy^ 



"p<JpUi n u n . 



(8) 



The operators {A, 7^,(7} satisfy the SU(2) quasispin algebra [21 
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A, A Y 



C/(2j 



C, A ] ] = 2A ] 



(9) 



A form similar to (|7]) was introduced by Lipkin et al. in their original 
paper fT3}j . The usual Lipkin model is obtained setting Ai = 0, since this term 
essentially renormalizes the single particle energy e. In our case it corresponds 
to the particular case Z = N = (2j + l)/2, which implies v p = v n = u v = u n , 
and x = K - However, we shall keep this term in order to be as close as 
possible to the realistic situation. 



3 Exact solutions 

To obtain the exact solutions of (0) we have performed a Holstein-Primakoff 
mapping |2l], |22| of this hamiltonian. It involves the substitution of pair of 
fermions by functions of the exact boson operators and b, which fulfill 
the exact commutation rule [b, b^} = 1. The relations between both set of 
operators are the following 



(10) 

The exact solutions are obtained using the boson basis 

V n b ] - 

The hamiltonian (0) has terms which change the number of bosons in 
two units, thus the exact wave function does not have a definite number of 
bosons. At the same time, states with odd and even number of bosons are 
not connected. To construct a decay scheme for double Fermi transitions 
the exact ground state of the even-even nuclei will be represented by the 
lowest energy state with an even number of bosons while the + states of the 
odd-odd nuclei are those with an odd number of bosons. Spurious states are 
avoided by limiting the number of bosons to < < 2j + 1. The physical 
states are written 
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2j+l 2j+l 

\Kven-even) = X/ ^n b \ n b) j |A dd-0(ta) = X/ C n °\n b ) . (12) 

"6 even r i b odd 

It must be clear that what we call "exact solutions" are the exact solutions 
of the hamiltonian (|7|), where the terms of the form ar at n , at^aip (scattering 
terms) have been neglected. These exact solutions exhaust any extension of 
the QRPA intended to diagonalize hamiltonian (0). However, the solutions of 
the complete eigenvalue problem of the hamiltonian (|l|) span a larger Hilbert 
space, including states orthogonal to those present in the actual "exact" 
solutions. The operator (3^ is sensitive to this truncation of the Hilbert space 
generated by solving the restricted problem defined by the use of hamiltonian 
(|7p. While at the QRPA level it has no effect, in any extensions beyond QRPA 
the transition amplitudes exhibit the missed components. Its effects on the 
Ikeda sum rule are described below. 



4 QRPA and RQRPA 

The QRPA hamiltonian Hqrpa can be obtained from Eq. (0) by taking the 
limit (2j + 1) — > oo. It is given by 

Hqrpa = (2e + Ax) 6+6 + A 2 { 6+6+ + bb}. (13) 

The QRPA states are generated with the one-phonon operator Oq RPA = 
XA< - YA acting over the correlated QRPA vacuum |0 >. The quasiboson 
approximation assumes the < 0\[A, A^]\0 >= 1, leading to the normalization 
condition X 2 — Y 2 = 1. The QRPA matrix is just a 2 x 2 one, with sub- 
matrices Aqkva — 2e + Ai and Bqkva — 2A 2 . The eigen-energy is Eqrpa = 
[(2e + Ai) 2 — 4A 2 .] 1 / 2 . It becomes a pure imaginary number if 2X 2 > 2e + Aj. 
It means that for this limit the zero-boson component of ground state ceased 
to be dominant. 

In the renormalized QRPA the structure of the ground state is included 
explicitly pTT | , in the form 



r -dtydt 

|0 >=Ne s \BCS> , S= , r „ . (14) 

1 1 ' 2<0p,^ t ]|0> 
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The RQRPA one-phonon state is given by 



OrqrpaI® >z 



Xtf - yA] I < Op, A*]\0 >^ 2 |0> 



(15) 



The condition O R q RP a\0 >= leads to the value c = y/X. After some 
algebra it is possible to show that < 0\[A, A*]\0 >= D = 1 - f^f [T2|, [TJ], 
and that 

r 2y 2 T 1 

b ^ 1 + ^ttJ ■ (16) 

The RQRPA submatrices are .At^q-ktm = 2e + XiD and B-jiqrva — 2A2-D. 
Given < -D < 1, the presence of D multiplying both Ai and A2 gives 
the needed reduction of the residual interaction to avoid the collapse of the 



QRPA equations fl3| . Due to this fact the RQRPA energy E R q RP a is always 
real. Its value must be obtained by solving simultaneously the non-linear 
equations for E R q RP a, X, y and D, which in the general case will include all 



possible values of J [R|. 



5 Results and discussion 

In the following we will present numerical results which correspond to the 
model space and parameters 



j = 9/2 , Z = 4 , N = 6 , e = IMeV . 



(17) 



In order to avoid dealing with small numbers, we now redifine the two pa- 
rameters 

k - k' = (2j + 1)k, X^X'= (2j + l)x- (18) 

We selected the values x' — or 0.5. The particle-particle strength k' is kept 
clS db variable. 



5.1 Excitation energies 

In Fig. 1 the excitation energies E\ — E , Eq RP a and E R q RP a are plotted 
against k' for x' — (upper figure) and x! — 0.5 (lower figure). In both 
cases the collapse of the QRPA is evident at k' ~ 1. The RQRPA excitation 
energy remains real but it decreases as compared with the exact one. For a 
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relatively large value of k' , («' ~ 2), the ground state and the first excited 
one tend to become degenerate. 
Fig. 1 

This figure strongly resemble Fig. 1 of ref. |13j and Fig. 2 of ref. [§], 



where the energy of the first excited state is plotted against the particle- 
particle strength parameter g pp . The curves for the QRPA and the RQRPA 
are quite similar to those shown here. The advantage of the present simple 
model is that we can compare them with the exact results, while in the general 
case with multiple single particle levels the exact solutions are unknown. 

5.2 Boson number 

The expectation value of boson number operator in the ground state measures 
the difference between the QRPA ( or RQRPA) ground states and the BCS 
vacuum. Its expected value should be half the number of quasiparticle and 
for the different cases it is written as 

2j+l 

<0\h\0> exact = £ \C° n \ 2 n, 

n even (19) 

< 0|n|0 > RPA = Y 2 , < 0|n|0 >rqrpa= y 2 ■ 

Fig. 2 shows the behavior of this average number for \' — and 0.5. Up 
to k' ~ 0.7 the three quantities are nearly indistiguishable. From there on, 



the QRPA overestimates this correlations, as pointed out long ago pi), and 
quickly collapses. The RQRPA does not collapse but it shows about twice 
the exact number of bosons. 
Fig. 2 

This behaviour will affect the Ikeda Sum Rule, as it will be shown below. 

5.3 (3 and (3(3 transition amplitudes 

The Fermi f}^ operators in the quasiparticle basis are 



(3 = y/Tj + 1 u p v n A^ + v p u n A - u p u n B^ + v p v n B 

/3+ = (/r)t, fit = [xjw 



(20) 

Neglecting the scattering terms B^,B, as it was done to obtain hamil- 
tonian ([?[), leads to the QRPA order of approximation. In this case the 
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amplitude of the Fermi transitions connecting the ground state |0) of the ini- 
tial even-even nuclei and the ground and excited states |0a > in the odd-odd 
nuclei are 



< Ox\P-\0 >e X act = 2 E Cl U p V n C: X +1 [(1 - 5^) {n + 1) 



n even 



VpU n C n _i 



2j 

1/2 



1,2 + 



2j+l ■ 



n 



< 0l|/3 |0 > RPA = y/2j +l (X u p v n + Y v p u n ) 

< 0l|/3 - |0 >RQRPA = a/ 2 J + 1 {X U P V n + y V p U n ) 

(21) 

Similar expressions hold for (3 + interchanging the w's and v's. There are 
several states in the odd-odd nuclei which can be connected to the ground 
state of the even-even one (five in our example). Only the state with the 
lowest energy, labeled with |0i >, is described by the QRPA and the RQRPA. 

In order to study the Pfov decay amplitudes M 2u in this simple model we 
made the approximation < /ina/ |/?"|0 A >^< in j tia /|/3~|0 A >=< 0\\f3 + \0 > 
[0, §]. In this way, we use the expression 

m 2 „ = £ < "^10 >< > ( 22) 

A ^A + A 

We selected A = 0.5 MeV. The sum runs over all the odd-odd states, 
which are only one in the QRPA and RQRPA. 

Fig. 3 shows M.2 V as a function of k' for x' — and 0.5. Their behaviour is 
very similar of that encountered in realistic calculations || ^ ^ [L3|, including 
its cancellation near the collapse in the QRPA description. The success 
of the RQRPA in extending this curve far beyond the value of k' where 
the QRPA collapsed is clearly seen. By the other side, compared with the 
exact values of M 2u , the RQRPA performance is poor. The overestimation 
of the ground state correlation causes a premature cancellation of the Pfov 
transition amplitude, as compared with the exact result. 

Fig. 3 

5.4 Ikeda Sum Rule 

There is an additional point which must be mentioned. The Ikeda Sum Rule 
S~-S + = N-Z, S ± = Y,\< OaI/^IO > | 2 (23) 

A 
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must be fulfilled in any model where the Hilbert space for the odd-odd nuclei 
includes all the states which can be connected to the ground state via the 
beta decay. 
Fig. 4 

The behaviour of S~ — S + against k' is exhibited in Fig. 4. As it is well 
known, in the QRPA this sum rule is always fulfilled, irrespective of how 
large the ground state correlations are. It is seen in the figure as a dotted 
straight line, interrupted where the QRPA collapses. The RQRPA strongly 
violates the sum rule by nearly 50 % at large values of k. The exact solution 
of (0) also violates the sum rule but in a smaller amount. The origin of 
this problem can be attributed to the adopted structure of both the hamilto- 
nian and the transition operators. As it was mentioned above neglecting the 
scattering terms in the operators H, (3~ and j3 + provides a plausible expla- 
nation of this failure. At the QRPA level these terms play no role, but when 
ground states correlations are explicitly taken into account they become rel- 
evant and cannot be neglected, neither in the transition operators nor in the 
Hamiltonian. 



6 Conclusions 

We have presented a solvable model which is holomorphic to the Lipkin 
model and used it to study the QRPA and RQRPA approaches. The case 
of double beta decay transitions of the Fermi type has been discussed. We 
have found that the excitation energies and double beta decay amplitudes, 
for these transitions, have the same qualitative behaviour found in realistic 
calculations. It is shown that, as expected, the RQRPA does not collapse as 
the QRPA description does when renormalized particle-particle correlations 
are included. However, the apparent stability of the RQRPA is hindered by 
the fact that it strongly overestimates the effect of ground state correlations. 
This tendency is reflected by the premature change in the sign of M 2u , as 
compared to the exact solution, and by the violation of the Ikeda Sum Rule. 
In our opinion it implies that if one wants to include explicitly the so-called 



ground-state-correlations in the fashion of |IT[ one should also exceed the 
two-quasiparticle order. The usually neglected one-quasiparticle terms are 
included in the RQRPA when the exact commutation relations for and A 
are used but they should be also included in the transition operators and in 



11 



the hamiltonian. In our exact model, the inclusion of these terms implies the 
use of the SO(5) algebra |22|, instead of the present SU(2) one. Details 
about the work in progress will be published elsewhere. 
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Figure Captions 

Fig. 1: Excitation energy vs. the particle-particle strength k'. The 
particle-hole strength is x' = (upper figure) and \' = 0.5. The exact 
solution is plotted with thin line, the QRPA with dots and the RQRPA with 
a thick line. 

Fig. 2: Number of bosons < 0|n|0 > in the ground state vs. k' . Conven- 
tions are the same as in Fig. 1. 

Fig. 3: PP 2 u transition amplitudes M 2u vs. k' . Conventions are the same 
as in Fig. 1. 

Fig. 4: Ikeda sum rule vs. k,'. Conventions are the same as in Fig. 1., for 
X 1 = 0.5 only. 
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Figure 1 
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